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ABSTRACT 

Type-IIn supernovae (SNe), which are characterized by strong interaction of their ejecta with the 
surrounding circumstellar matter (CSM), provide a unique opportunity to study the mass- loss history 
of massive stars shortly before their explosive death. We present the discovery and follow-up observa- 
tions of a Type Iln SN, PTF09uj, detected by the Palomar Transient Factory (PTF). Serendipitous 
observations by GALEX at ultraviolet (UV) wavelengths detected the rise of the SN light curve prior 
to the PTF discovery. The UV light curve of the SN rose fast, with a time scale of a few days, to a 
UV absolute AB magnitude of about —19.5. Modeling our observations, we suggest that the fast rise 
of the UV light curve is due to the breakout of the SN shock through the dense CSM (n « 10 10 cm -3 ). 
Furthermore, we find that prior to the explosion the progenitor went through a phase of high mass-loss 
rate (~ 0.1 Mg yr _1 ) that lasted for a few years. The decay rate of this SN was fast relative to that 
of other SNe Iln. 

Subject headings: stars: mass-loss — supernovae: general — supernovae: individual (PTF 09uj) 



1. INTRODUCTION 

Early ultraviolet (UV) detection of supernovae (SNe) 
of all types in the shock-breakout phase holds great po- 
tential for probing the nature and properties of SN pro- 
genitors (e.g., Colgate 1974; Falk & Arnett 1977; Klein 
& Chevalier 1978; Matzner & McKee 1999; Waxman ct 
al. 2007; Katz, Budnik & Waxman 2010; Piro, Chang & 
Weinberg 2010; Rabinak & Waxman 2010; Nakar & Sari 
2010). Shock-breakout observations provide a measure 
of the progenitor radius and the structure of its outer 
layers (e.g., Soderberg et al. 2008; Gezari et al. 2008; 
Schawinski et al. 2008). 

Type Iln SNe are characterized by the presence of a 
blue continuum and narrow emission lines in their opti- 
cal spectra (e.g., Schlegel 1990; Filippenko 1997). These 
features are usually interpreted as the signatures of inter- 
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action of the SN ejecta with dense circumstellar matter 
(CSM), due to stellar mass loss prior to the explosion. 
SNe Iln are probably an inhomogeneous class of objects 
whose main properties are dictated by the presence of 
dense CSM rather than by the details of the explosion. 
Therefore, SNe showing evidence for strong interaction 
are unique probes of scenarios in which a stellar explosion 
follows a major mass-loss event, and they can be used 
to study the mass-loss rate from progenitor stars (e.g., 
Chugai & Danziger 1994; Chugai et al. 1995; Fransson 
et al. 2002; Gal- Yam et al. 2007; Ofek et al. 2007; Smith 
et al. 2007; Smith & McCray 2007; Gal- Yam & Leonard 
2009; Dessart et al. 2009). 

Here we present the discovery of a Type Iln SN, 
PTF 09uj, which was serendipitously observed in the UV 
by GALEX (Martin et al. 2005) shortly after the explo- 
sion. 

2. OBSERVATIONS AND DATA REDUCTION 

SN PTF09uj was discovered on 2009 June 23.30 
(UTC dates are used throughout this paper) by the 
Palomar Transient Factor^ (PTF; Law et al. 2009; 
Rau et al. 2009) conducted with the Oschin 48- 
inch Schmidt telescope (P48) at Palomar Observatory. 
The SN is associated with an r = 17.3 mag galaxy, 
SDSSl3j142010.86+533341.9 (Fig. HJ — a late- type disk 
galaxy at a redshift of z = 0.0650 ± 0.0001 (distance 

d = 292Mp<E!D- 

2.1. Photometry 

Follow-up photometry was obtained using the auto- 
mated Palomar 60-inch telescope (Cenko et al. 2006; 

13 http:/ /www. astro. caltech.edu/ptf/. 

14 Sloan Digital Sky Survey; York et al. (2000). 

15 We assume WMAP-5 cosmological parameters (Komatsu et 
al. 2008) and insignificant peculiar velocity. 
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TABLE 1 
Observations of PTF 09uj 



Fig. 1.— Images of the field of PTF09uj. The SN was de- 
tected by GAL EX in the NUV band on June 22 (upper left), 
and by PTF in the R band on June 23 (upper right). The 
bottom-right panel shows a reference ij-band image prepared from 
PTF images obtained before the explosion. The numbered stars 
(1—6) mark the SDSS photometric reference stars. The bottom- 
left panel shows the discovery of PTF 09uj using image sub- 
traction of PTF data. A 30"-radius circle marks the position 
of the SN in all panels. The position of the SN is a(J2000) 
= 14 h 20 m 11.15 s , 5(J2000) = +53°33'41.0", which is 2.7" from 
the center of SDSS J142010. 86+533341. 9 at a position angle of 

110 deg. 



Table [T]). As noted above, this field was observed by 
GALEX on several occasions (Table [1]). 

Photometry of PTF09uj in the PTF and P60 im- 
ages was obtained by the common point-spread function 
image-subtraction method (Gal- Yam ct al. 2004; Gal- 
Yam et al. 2008). Estimates of uncertainties were ob- 
tained from the scatter in the magnitudes of artificial 
sources. The photometry was calibrated using r-band 
magnitudes of six SDSS stars (Fig. [T]). Calibration er- 
rors were added in quadrature to the image-subtraction 
errors (Table [lj. 

The GALEX photometry was carried out by perform- 
ing aperture photometry with a 10" radius around the SN 
host galaxy, and subtracting its light as measured in the 
reference image. The reference image was constructed 
by combining the four GALEX images of this field taken 
prior to the SN explosion, between 2009 May 4 and May 
14. Since GALEX uses photon-counting detectors (i.e., 
individual photons are time tagged; Martin et al. 2005), 
we had the opportunity to look for flux variations on rel- 
atively short timescales. In particular, we examined the 
earliest image in which the SN was detected, which was 
taken on June 22 and had an exposure time of 1364 sec- 
onds. We extracted the time-tags of the 4597 photons 
found within 10" of the SN and binned these photons 
on timescales from 3 to 1000s. We found no significant 
variations in flux as a function of time. We note that 
about 40% of these photons originate from the SN and 
the rest are due to the host galaxy. The GALEX NUV 
and P48 i?-band light curves of PTF 09uj are presented 
in Figure [2] 
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Note. — /„ is calculated at 2316 A, 4718 A, 6184 A, and 7499 A 
for the NUV, g, PTF R/r, and i bands, respectively. Magnitude 
uncertainties include (in quadrature) absolute calibration errors of 
0.099 mag for the PTF i?-band measurements, and 0.071, 0.037, 
and 0.029 mag for the P60 g, r, and i bands, respectively. An 
aperture correction of 0.12 mag was applied to the GALEX- NUV 
magnitudes (Morrissey et al. 2007). 

a The last P48 non-detection before the discovery. 
b All of the P48 observations were conducted using the Mould R- 
band filter. Photometry was measured in the combined images of 
the same field taken each night (usually two). 

c The magnitude from Lick observatory is based on synthetic pho- 
tometry of the spectrum using the code described by Poznanski et 
al. (2002). 
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Fig. 2.— The light curve of PTF09uj from the P48 (red circles), 
P60 (blue diamonds), synthetic photometry (magenta star), and 
GALEX NUV observations (black squares). The tip of the black 
triangle marks the magnitude lower limit from GALEX. The curves 
show the scaled and smoothed light curves of two other SNe Iln 
(SN1998S, Fassia et al. 2000; SN 2005gl, Gal- Yam et al. 2007). 



We observed PTF 09uj with an exposure time of 1800 s 
with the Kast double spectrograph (Miller & Stone 1993) 
mounted at the Cassegrain focus of the Shane 3-m tele- 
scope at Lick Observatory; the 5500 A dichroic was em- 
ployed. On the blue arm, we used the 600 lines mm -1 
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Fig. 3. — The spectrum of PTF 09uj (black line; upper spectrum) 
and its host galaxy (red line; lower spectrum). For comparison, the 
scaled spectrum of the Type Iln SN 1998S (Fassia et al. 2001) two 
days prior to maximum light is shown just below the spectrum of 
PTF09uj (gray line). The inset shows the Ha line in detail. 

grism blazed at 4310 A to provide spectral coverage of 
3500-5550A and a dispersion of 1.02 Apixel -1 , while on 
the red arm we used the 300 lines mm -1 grating blazed 
at 7500 A for a wavelength range of 5400-10000 A and a 
dispersion of 4.60 Apixel -1 . 

The Lick spectrum was reduced using standard rou- 
tines in IRAIO (details provided by Cenko et al. 2008). 
Flux calibration was performed relative to the standard 
stars BD +28° 4211 (blue side) and BD +26° 2606 (red 
side). The Lick spectrum of PTF09uj is shown in Fig- 
ure [H Given the lack of Na I absorption lines in the SN 
spectrum and the low Galactic extinction toward this 
SN (E B ~v = 0.011 mag; Schlegel et al. 1998), we do not 
correct for extinction. 

3. INTERPRETATION 

The spectrum of PTF09uj, which was taken around 
peak light, shows a blue continuum, with narrow Ha 
emission and no prominent, broad absorption features. 
This "narrow" line is actually broader than the Ha line 
from the host galaxy and shows a hint of a P-Cygni pro- 
file (Fig. [3] inset). The spectrum also exhibits a narrow 
He I emission line (measured rest wavelength 5884 A, cor- 
responding to He I A5876) which is not present in the 
host-galaxy spectrum. These observations suggest that 
PTF 09uj is a SN Iln enshrouded with a dense CSM. 

However, the e-folding decline rate of the SN flux is 
about 10 days. This is faster than the steepest declining 
SNe Iln previously known, such as SN 1998S (Fassia et al. 
2000), SN2005gl (Gal- Yam et al. 2007), and SN2005ip 
(e.g., Smith et al. 2009); compare with the light curves 
of the first two in Figure [5] 

Another possible difference between PTF09uj and 
SN 1998S is the spectra. While the SNe spectra shown 
in Figure |3] were both taken around maximum light, 
PTF 09uj evolved faster and therefore these spectra prob- 

16 IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association for Research in 
Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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Fig. 4. — A sketch of our model for PTF09uj. Fast ejecta from 
the SN explosion are interacting with a CSM having a wind-like 
density profil e. Here i> w is the wind speed. See definitions of vari- 
ables in 93.11 

ably do not correspond to the same epoch after explosion. 
In order to compare the spectra of the SNe taken at the 
same epoch after explosion it is probably more adequate 
to inspect the spectrum of SN 1998S taken 13 days prior 
to maximum light (Fassia et al. 2001). This earlier spec- 
trum of SN 1998S is different, with broader emission lines 
and strong "Wolf-Rayet"-like features (e.g., C III, N III). 

3.1. Shock Breakout in a Stellar Wind 

The fast rise in UV light and the high peak luminos- 
ity (vf u r; 3 x 10 43 ergs _1 ) motivates us to consider a 
model of a shock breakout which takes place within a 
dense, optically thick, stellar wind (see also Falk & Ar- 
nett 1977; Waxman et al. 2007). The blue continuum in 
the visible-light spectrum of PTF 09uj suggests that the 
emission is optically thick. Fitting a black-body curve to 
the NUV, g, r, and z-band photometry, obtained on June 
25-26, we find a best-fit temperature of « 1.7 x 10 4 K 
(with r.m.s. of 0.13 magnitudes). We note that if line 
blanketing is effecting the spectrum, the true effective 
temperature could be even higher. 

In the framework of the model considered here (see 
sketch in Fig. 2]), the rising UV emission is due to a 
shock breakout within an optically thick wind. Some or 
most of the visible-light emission at later times is caused 
by diffusion of the shock-deposited energy. 

We now calculate the properties of the shock and ejecta 
that are needed in order to explain the observations. In 
this calculation we use the observed peak luminosity and 
the rise time of the SN to calculate various parameters 
(i.e., mass, velocity, temperature). 

The thickness of a radiation- mediated shock, r s , in 
units of the Thomson optical depth (e.g., Weaver 1976) 
is given by 

r s « c/v s . (1) 

Here c is the speed of light and v s is the upstream ejecta 
(and shock) speed. A radiation-mediated shock "breaks 
down" or "breaks out" (i.e., radiation escapes ahead of 
the shock) when photons diffuse ahead of the shock faster 
than the shock propagates. For a wind-density profile 
p(r) — Kr~ 2 where if is a normalization constant and 
r is distance from the center, the photon diffusion time 
from r to 2r is t^is = npr 2 /c = kK/c, independent of 7', 
while the shock propagation time is r/v s , growing with 
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r. Thus, the shock breaks down when it reaches rb rca k = 
kKv s /c. At this point, photons would diffuse and escape 
the wind on a time scale equal to idiff multiplied by a log 
correction factor ln(c/w s ). 

An alternative definition of the breakout radius is the 
radius at which the thickness of the shock is comparable 
to the scale height of the density variation. Since for a 
wind profile this scale is r, with a corresponding optical 
depth r = upr — nK/r, this happens at r s = nK/r = 
c/v s or rbrcak = kKv s /c, the same as obtained above. 
For small radii r, the shock expansion time scale (t exp = 
r/v) is shorter than tdiff, and breakout takes place when 
idiff ~ ^exp- We further note that v s corresponds to the 



occurs is 



velocity of the faster parts of the shock. In £13.21 we will 
derive the ratio of this velocity to the bulk velocity of 
the ejecta («b)- 

Assuming a radiation-dominated shock, the down- 
stream temperature of the shock can be estimated by 
comparing the radiation energy density (assuming black- 
body radiation) with the kinetic energy per unit volume: 



aT 4 



7 



7c 
27k' 



(2) 



where a is the radiation constant, t is the time since 
the explosion, k is the opacity, and p s is the density of 
the stellar wind at the shock breakout. Note that the 
right-hand side of Equation[2]is derived assuming a wind- 
density profile p(r) oc r~ 2 along with Equation [1] and 
r = J Kp(r)dr = Kpr s . The factor of 7 arises from the 

shock's compression ratio (7+ l)/(7 — 1) = 7, with an 
adiabatic index of 7 = 4/3. Here r s is the radius at which 
the shock breakout takes place. 

In this framework, the kinetic energy of the explosion 
is converted into radiation at a rate of L 2itr 2 v s p s v 2 . 
Here we assume that the shock is optically thick and that 
the emission spectrum is roughly represented by a black- 
body spectrum. The assumption of black-body emission 
is justified, since for the shock velocity we infer, v/c ~ 
0.03 (see Eq.|4]below), there is no departure from thermal 
equilibrium. Such departure is expected only for higher 
velocity shocks, v/c i> 0.1 (Weaver 1976; Katz et al. 
2010). 

Based on Equation [TJ and assuming we are observing 
at wavelength A in the Rayleigh- Jeans tail, 



L x » —L 



BB 



2 2 2k B T 

47r r s Vs— 



(3) 



where L^ B is the black-body total specific luminosity 
and kg is the Boltzmann constant. 

By solving Equations [5] and [3J along with r = v s t, we 
find that 



v ~ 1 3 x 10 4 L 1/3 A 4/3 «- 1/12 r 7/12 km s" 1 

U s ~ 1.0 X 1U 1^6038^6200^0.34 c 7 Km h ' 



(4) 



where £,\,6c38 * s ^ ne specific luminosity in units of 6 x 
10 38 ergs -1 A -1 (equivalent to the luminosity measured 
on 2009 June 23, close to the first GALEX detection), 
A.6200 is the wavelength at which the luminosity is mea- 
sured in units of 6200 A, K0.34 is the opacity in units of 
0.34 cm 2 g _1 (assuming the opacity of completely ionized 
solar composition) , and tj is the time between the explo- 
sion and the measurement of L\ 1 in units of 7 days. The 
distance from the star to the shell in which the breakout 



1, t ~ 8 9 x If) 14 / 1/3 A 4/3 k 1/12 / 5/12 , 
U S L ~ O.Z X 1U -^x 6c38 A 6200 K 0.34 L l ' 



(5) 



This radius is larger than the size of a typical red su- 
pergiant, indicating that our assumption that the shock 
occurs "outside" the star, in a wind-density profile, is 
justified. The density is given by 



xnn 10r-2/3 ,-8/3 -7/6.1/6 -3 
x iu ^A,6c38 A 6200 K 0.34 L 7 cm > 



(6) 



where m p is the proton mass, and we also have 

T»9.1 x 10 4 k~3/V 1/4 K. (7) 

Here, T is the effective temperature at the time in which 
we measured L\. 

In addition, we can estimate the rise time of the SN 
light curve. In the case that the scale height of the ma- 
terial lying ahead of the shock is negligible compared 
with the radius at which the shock occurs (r s ; i.e., shock 
takes place at the edge of the star), the rise time is ~ r s /c 
(up to a log correction factor mentioned earlier). How- 
ever, in the scenario discussed here there is a significant 
amount of material ahead of the shock and the rise time 
is i r i sc r s /v s = 7t? day. This rise time is therefore 
slower than the rise time seen in some of the known shock 
breakout events (e.g., Gezari et al. 2008). In this sce- 
nario, after maximum light, the light curve is expected 
to decay exponentially on the diffusion time scale (idiff; 
Falk & Arnett 1977), which in our case is about a week. 
This is similar to the observed decay rate. 

The total mass in the fast ejecta within radius r s is 
J r " 4irr 2 p(r)dr . Here r* is the star radius. Assuming the 

wind density is p(r) oc r~ 2 and r/r* 3> 1, we find 

A S^noy-i/S \4/3 -11/12.17/12 / s 

Airpri w 0.3L a ( 6o38 A 6200 k u ' t 7 M . (8) 



M 



Note that the total mass in the ejecta may be larger (see 
£|3.2p . and therefore this equation provides only a lower 
limit to the mass. Next, the mass loss rate from the 
progenitor prior to the explosion is 



M 



Mv. 



Ml, 100 



0.1k 



0.34*7^,100 M © y~ 



(9) 



where v^ioo is the progenitor wind velocity in units of 
100 km s -1 . This mass-loss rate is required to persist for 
Ts/vw ~ 10 yr prior to the explosion. 

Another important property is the kinetic energy in 
the faster parts of the shock, which is roughly given by 



c Y in50r \4 -3/4.1/4 

O X 1U 1jA : 6c38 a 6200 k 0.34 l 7 



erg. (10) 

As before the factor of 7 arises from the shock's compres- 
sion ratio. 

We note that if we relax the Rayleigh- Jeans approx- 
imation in Equation G2 and solve these equations using 
the full Planck formula, the solution does not change con- 
siderabl"vF1. Therefore, for simplicity, we choose to show 
here the approximate solution. 

17 The values of v s , r s , t r i se , and M are changed by a factor of 
1.04; n is changed by a factor of 0.92; E by 1.14; and T and M 
remain the same. 
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Based on the temperature that we derive in Equation[?l 
the expected NUV luminosity is 2.1 x 10 40 ergs -1 A -1 
(given by the Planck function multiplied by v s /c). This 
is about 1.5 times larger than observed. However, the UV 
emission may be affected by metal line blanketing, and 
is therefore less reliable than the visible-light luminosity. 

Furthermore, we note that our assumption of a 
radiation-mediated shock is justified since in the 
case studied here the ratio of radiation energy den- 
sity to plasma thermal energy density is very large, 
aT 4 /{nk B T) w 10 6 . 

To summarize, based on our crude model we use the 
visible-light luminosity to derive the shock properties 
(e.g., r s , v s , T, E, M). The calculated t rise is roughly 
consistent with the observed rise time of the UV light 
curve. Moreover, this model naturally explains the high 
observed luminosity in the NUV band. 

3.2. The Bulk Velocity 

The velocity v s we used so far corresponds to the faster 
parts of the ejecta. However, this does not necessarily 
represent most of the energy in the ejecta. As a sanity 
check, here we estimate the bulk velocity of the ejecta. 

The shock that accelerates the ejecta gives more energy 
to slower shells. The energy as a function of velocity is 
given by 

(11) 



E(v)=E h (— 



where E b is the bulk kinetic energy of the ejecta, v\> is the 
bulk velocity, and x = 5(l + 3n/5)/n, where n = 3/2 (3 
x = 19/3 (14/3), for convective (radiative) envelope: 
(Matzner & McKee 1999). A shell of velocity v is de- 
celerated when the energy in the shocked wind, E w , 
equals the energy in the decelerated shell, E(v). Here 
E w = 4nKrv 2 and E(v) = M(v)v 2 . Note that this in- 
cludes the internal energy in the shocked wind (assuming 
the internal energy roughly equals the kinetic energy). 
This gives the deceleration radius 



r(v) 



E h 



4irKv 2 V Vh 



The wind optical depth at this radius is npr, or 
AtthK 2 vI I v \ 



t(v) 



E h 



\v y 



(12) 



(13) 



As long as r ^ c/v s , the fast shell is decelerated and 
overtaken by the slower, more energetic, shells behind it. 
Breakout occurs when t(v) ~ c/v s , which together with 
Equation [TBI and assuming Eh = M b v 2 /2 gives 



Hi - ( — 

Vb V« s 



l/(2+x) 



8ttkK 2 J 



(14) 



where Mb is the total mass in the ejecta. Since l/(x + 
2) <C 1, the inferred value of v s /v b is insensitive to the 
exact values of v s and K, and to the unknown value 
of M b . Assuming M b = 1M Q , 

K = M/(4nv w ) w 5 x 10 16 



g cm 



v s = 10 4 kms 1 , and 
1 (corresponding to 



18 In the notation of Matzner & McKee (1999), x = (n + 
l)/(n/3i) - 2, where /3i = 1/5. 



M = 0.1M©yr -1 and v w = 100 kms -1 ), we get v s /v b w 
1.1 (1.2) for convective (radiative) envelopes. 

This analysis suggests that Equations U and [TO] are 
reasonable approximations of the bulk properties of the 
ejecta. Specifically, if the total mass of the ejecta is an 
order of magnitude larger than that given by Equation[Sl 
then the total kinetic energy of the SN will exceed 5 x 
10 51 erg, which is unlikely (at least for the garden variety 
of SNe). 

4. DISCUSSION 

We present the discovery of PTF09uj, which was 
serendipitously observed by GALEX at early times after 
the explosion. The spectrum of the SN and the bright 
UV signal suggest that this was a SN Iln powered by the 
diffusion of the shock energy and interaction of the ejecta 
with a dense CSM (n « 10 10 cm -3 ). This interpretation 
is consistent with both the fast rise of the UV light curve 
and the UV luminosity. The observed fast rise cannot be 
easily explained unless the progenitor is embedded in an 
optically thick wind. Moreover, a shock breakout from 
a stellar photosphere cannot generate such a bright UV 
signal (e.g., Rabinak & Waxman 2010). 

Based on simple modeling, we suggest that prior to the 
explosion the progenitor went through a phase of high 
mass- loss rate, with M 0.1 M Q yr -1 . The radius of the 
radiating region and the fast decay of the SN are sugges- 
tive of an episodic high mass-loss rate with a duration 
of about several years prior to the explosion. Our model 
suggests that the total mass of the ejecta is relatively 
low, roughly 1 M Q . We stress that this is an order-of- 
magnitude estimate; the total ejected mass could be as 
high as a few solar masses, but probably not on the high 
end (^ 10 M ) of the ejected mass expected in typical 
SNe II. 

The low ejecta mass, if true, may be due to one of the 
following: (i) most of the mass of the progenitor was not 
expelled by the SN explosion and it is in a form of a com- 
pact remnant; or (ii) most of the mass of the progenitor 
was expelled (e.g., wind) prior to the explosion. This is 
in contrast to more energetic SN explosions whose lumi- 
nosity is powered by interaction with a dense CSM (e.g., 
Benetti et al. 2006; Ofek et al. 2007; Smith et al. 2007) 

PTF is a wide and shallow survey in which about 
3000 deg 2 are actively surveyed at a given time down 
to a limiting magnitude of ~ 21. GALEX sensitivity 
allows high signal-to-noise ratio detections of SN shock- 
breakout flashes at a redshift range which is similar to 
that probed by PTF. Specifically, we estimate that PTF 
should find several SNe each year for which GALEX early 
observations will be available. This estimate is based on 
the number of PTF SN discoveries, during 2009, which 
had GALEX serendipitous observations between to 30 
days prior to the PTF discovery. 
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